Tunneling conductance and local density of states in time-reversal symmetry breaking 
superconductors under the influence of an external magnetic field 
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We consider different effects that arise when time-reversal symmetry breaking superconductors are subjected to 
an external magnetic field, thus rendering the superconductor to be in the mixed state. We focus in particular 
on two time-reversal symmetry breaking order parameters which are believed to be realized in actual materials: 
p + lp'-wave and d + is- or d + id'-wave. The first order parameter is relevant for S^RuOzt, while the latter 
order parameters have been suggested to exist near surfaces in some of the high-T c cuprates. We investigate the 
interplay between surface states and vortex states in the presence of an external magnetic field and their influence 
on both the tunneling conductance and the local density of states. Our findings may be helpful to experimentally 
identify the symmetry of unconventional time-reversal symmetry breaking superconducting states. 
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I. INTRODUCTION 



Recently, considerable attention has been devoted to the 
chiral superconducting phase which is believed to be real- 
ized in the p-wave triplet superconductor 1 Sr2Ru04. The chi- 
ral state of a p-wave superconductor corresponds to a non- 
zero projection l z — ±1 of the Cooper pairs angular mo- 
mentum 1 along the z axis, and thus breaks time-reversal 
symmetry. The spatially homogeneous triplet order param- 
eter A = Ao(d ■ a)ia y is described by the vector 1 d(p) = 
(0, 0,p x + ixPy)' which depends on the direction of electron 
momentum p. Here Ao is the bulk value of the order pa- 
rameter, a = (a x , a y) a z ) is the vector of Pauli matrices of 
conventional spin operators, and x = ±1 corresponds to the 
two possible values of chirality. Also, chiral superconduct- 
ing states can be associated with an admixture of two order 
parameters corresponding to different irreducible representa- 
tions of a crystal point group. Naturally different order param- 
eter components can coexist in the vicinity of interfaces be- 
tween superconductors and surfaces due to broken symmetry 
of the crystal group 2,3 ' 4 ' 5 . Among other possibilities of sub- 
dominant order parameter symmetry 2 there are states which 
break time reversal symmetry 3 - 4 - 5 . The coexistence of order 
parameters shows up in the local density of states 6 - 7,8,9 as well 
as in the generation of spontaneous currents flowing along the 
surfaces in time reversal symmetry breaking cases 7 -^. 

Time-reversal symmetry breaking order parameters have 
been proposed to exist near surfaces 1 - and within vortex 
coresi 1 in high-T c superconductors. This proposal stems from 
the observation of a split zero-bias conductance peak in the 
absence of any applied magnetic field. In this case, it has 
been suggested that the relevant order parameter is either 
d + is- or d + id'- wave. The gap may then be written as 
A = A g(6 p ) + iA s or A = A g(8 p ) + iA d gx(9 p ), respec- 
tively, where Ao is an amplitude of the main component and 
the admixture of another pairing symmetry is denoted by the 
amplitudes A s and A<j. Here, 8 p is a polar angle in momen- 



tum space p = p(cos 8 p , sin 8 p ), g(6 p ) — cos(26* p + a) and 
gi(9p) — sm(28 p + a) where a/2 is an angle measuring the 
misorientation of crystalline symmetry axes and coordinate 
axes. One obtains c? a .2_ J/ 2-wave symmetry of the main order 
parameter component for a = and d^-wave pairing for 
a = 7r/2. While the experimental data so far clearly indicates 
an order parameter which breaks time-reversal symmetry, the 
question of whether the symmetry is d + is- or d + trf'-wave 
remains unsolved. Clearly, experimental signatures that may 
distinguish these two types of pairings would be highly desir- 
able. 

One of the important features of unconventional supercon- 
ductors is the possibility for the existence of surface Andreev 
bound state a 12 ' 13 ' 14 . They occur in the vicinity of the scatter- 
ing interface between a superconductor and an insulator, if the 
incident and reflected quasiparticles (QP) with different mo- 
mentum directions see different phases of the order parameter. 
The consequence of the Andreev bound states formation is an 
increase of the local density of states at the surface resulting in 
zero-bias conductance peak anomaly observed 18 in tunneling 
spectroscopy of high-T c cuprates with d-wave symmetry of 
superconducting pairing as well as in 19 the p-wave triplet su- 
perconductor Sr2RuO_4. Also, the Andreev bound states deter- 
mine the anomalous low-temperature behaviour of the Lon- 
don penetration length 15 and the Josephson critical current in 
d-wave 16 and chiral superconductors^ 7 -. 

Under the influence of an applied magnetic field, screen- 
ing currents and vortices may be generated in a supercon- 
ductor. As a result, the spectrum of surface states acquires 
a Doppler shift, leading to a splitting of the zero-bias conduc- 
tance peak 10 . Abrikosov vortices located near a superconduct- 
ing surface generate an essentially inhomogeneous superfluid 
velocity field, which leads to a non-trivial electronic structure 
of the surface bound state s 20 ' 21 ' 22 . Also, it was recently^ 2 - pro- 
posed that the same Doppler shift effect should lead to a chi- 
rality selective influence of the magnetic field on the surface 
states in a p-wave chiral superconductor with broken time- 
reversal invariance. The quasiparticle density of states (DOS) 
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near a flat surface was shown to depend on the orientation of 
magnetic field with respect to the chirality as well as on the 
vorticity in case when the Abrikosov vortex is pinned near 
the surface of superconductor. Additionally, in superconduc- 
tors featuring gap nodes, such as in the case in pure d x 2_ y 2- 
symmetric superconducting cuprates, a vanishing pair poten- 
tial in nodal directions results in important ramifications for 
the physics of the system 1 1,23 ' 24,2S - 26 . 

To understand the effect of an externally applied mag- 
netic field on the surface DOS, let us consider a spectrum 
of Andreev bound states near a flat surface of a time-reversal 
symmetry breaking superconductor occupying the half-space 
x > 0. Below, we focus on the p + lp-wave, d xy + is- and 
d X y + id x 2_ y 2 -wave cases for concreteness. We consider a 
model situation assuming spatially homogeneous gap func- 
tion, having the following form in momentum space: 

A = A e lxe " (1) 

for p + ip-wave, 

A = A sin(26> p ) + iA s (2) 
for d xy + is- wave and 

A = A sin(26» p ) + iA d cos(20 p ) (3) 

for d X y + id x 2_ y 2-wave superconductors. 

Assuming that the QP are specularly reflected at the surface 
of the superconductor within a Doppler shift approach 2 ^, the 
spectrum of the surface states can be expressed as follow o 17 ' 29 : 
£a = £a0 + tdi where e a0 is a position of energy level in zero 
magnetic field and ed = Kkpv s is the Doppler shift energy 
which is determined by a local field of superfluid velocity. 
The superfluid velocity v s near a surface has only a tangential 
component, directed along the y axis v s = (0,v sy ,0), and 
can be related to the density of supercurrent flowing along 
the surface j s = env s , where e is electron charge and n is a 
concentration of Cooper pairs. The magnetic field is screened 
in a superconductor at the London length A as follows B = 
He~ x / X , where H is a value of magnetic field outside the 
superconductor. Therefore, the superfluid velocity is v sy = 
-(2e/mc)XH. 

If the magnetic field is absent, the spectrum of surface states 
is given byi 2 ^ 

e 0a = X^ak y / k F (4) 

for a chiral p-wave, 

e a0 = A s sgn(fc y ) (5) 

for a d xy + is-wave and 

e a0 = A d sgn(k y ) cos(26» p ) (6) 

for a d xy + %d x 2_ y 2 - wave. Here k y is a projection of QP mo- 
mentum along the surface. The above spectra may be formally 
obtained by solving-^ 

A(fl p ) = A(tt - Op) 

£aO -V|A(0 p )| 2 -^ o e aQ + W\A(Tr~8 p )\2-ei 

(7) 



In the d x 2_ y 2 + is wave case, one finds that e a o = 
■J Ag cos 2 {20 p) + A 2 , from which one infers that there are 
no subgap surface states. This is qualitatively different from 
the d xy + is wave case. The interesting effects occur in the 
latter case, so we focus on the d xy + is{d) wave symmetry in 
the following, corresponding to a = n/2. 

The transformation of these spectra due to the Doppler shift 
effect is shown in FigQ] To be definite we assume that A s > 
0, Arf > and % = 1. Considering the DOS at Fermi level, 
v = | de a /dk y | "Jq, in a chiral p-wave superconductor one can 
see that its dependence on the magnetic field is mono tonic: it 
either increases or decreases for different field directions (see 
FigOJ) as discussed in Ref. |22| 

Another behaviour of the DOS occurs in the case of a d + is 
wave superconductor. From Fig. QJ it follows that for a cer- 
tain field direction there are no states at the Fermi level e = 
(red dashed lines in FigQJi). For the opposite field direction 
(blue dash-dotted lines in FigJTJ)), intersections of spectral 
branches with the Fermi level appear when the superfluid ve- 
locity is large enough \v sy \ > A s /pp so that the value of 
momentum projection at the intersection point is smaller than 
the Fermi momentum \k*\ < kp. Thus, one can expect that 
the DOS at the Fermi level should be zero when H < H*, 
where H * is the magnetic field value providing the condition 
\v sy \ = \A s \/pf to be fulfilled. 

On the contrary, in the d + id wave case the DOS at the 
Fermi level is non-zero even in the absence of a magnetic 
field. As can be seen from Fig. [TJ; (black solid lines) the 
spectral branches intersect the level e = at k* = ±kp/\/2. 
The transformation of the spectrum due to the magnetic field 
of different directions is shown in Fig. [TJ; by red dashed lines 
(H > 0) and by blue dash-dotted lines (H < 0). Then, it can 
be easily seen that for H > the coordinates of the intersec- 
tion points k y shift towards ±kp and for a certain value of the 
magnetic field H > H* the DOS at the Fermi level e = 
disappears. 

In the presence of an Abrikosov vortex near the surface of 
chiral superconductor a non-trivial structure of the local den- 
sity of states distribution appears which depends on the vortex 
orientation 22 . Along with the Doppler shift effect 22 , an impor- 
tant modification of the quasiparticle spectrum and the DOS 
can be obtained due to the overlapping of the surface states 
and the low-energy QP states localized within the vortex core, 
found in the pioneering work by Caroli- de Gennes and Ma- 
tricon (CdGM) 3 1 It was shown that QP states with energy 
lower than the bulk superconducting gap value A are localized 
within the vortex core at the characteristic scale of the order 
of coherence length £ and have a discrete spectrum e v (/i) as 
a function of the quantized (half-integer) angular momentum 
/i. At small energies |er| <gC A the spectrum for a vortex with 
vorticity M is given by 

e w (/i) w -M(j,uj, (8) 

where kp = pp/H and uj ~ Ao/{kp£). For the most of 
superconducting materials, including Sr2Ru04, the interlevel 
spacing ui is much less than the superconducting gap A since 
{kp£,) ^> 1. Therefore, the CdGM spectrum may be consid- 
ered continuous as a function of the impact parameter of the 
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FIG. 1 : (Color online) Plot of the surface states spectrum for (a) chi- 
ral p + ip-wave superconductor, (b) d + is- wave and (c) d + id- 
wave superconductors. Spectrum in zero magnetic field is shown by 
solid lines. Blue and red lines correspond to the spectrum transfor- 
mation due to magnetic field directed along and opposite to z axis 
correspondingly. 



quasiclassical trajectory b = —p,jkp and the direction of QP 
momentum ( 



Jp as in the following form 

£„(Mp)aMA(i/0. 



(9) 



In the case of a chiral p-wave superconductor, the spectrum of 
vortex core states differs from the CdGM result and is given 
by Eq. ([8]) with integer fi. For the d + is- and d + id- wave 
superconductors the quasiclassical spectrum of vortex core 
states is given by Eq.® with A(9 p ) = V /A^ 2 (6' p ) + and 
A(0p) = yj&lg 2 (6 p ) + A 2 ,g 2 (9 p ). The discrete spectrum 
is obtained by applying Bohr-Sommerfeld quantization rule 
to the canonical variables /j, = —kpb and 9 p 23 ' 31 . It should 
be noted that when the superconducting order parameter con- 
tains nodes, the quasiclassical expression (0 is invalid near 
the nodal directions since energy states near the vortex core 
are not truly localized, but rather "leak" out through the gap 
nodes 1 'i 23 ! 24 ' 25 ! 26 . This is not the case for us since we consider 
superconducting order parameters which are gapped over the 
entire Fermi surface. 

To study the interaction between vortex and surface states, 
let us consider an example of vortex positioned near a flat sur- 
face of chiral p + ip-wave superconductor. Comparing the en- 
ergies of surface e a dU and vortex e v © states one can see that 
for certain QP trajectories the condition of resonance e a = £ v 
is realized. Thus the spectrum transformation in such almost 
degenerate two-level system is given by a secular equation: 



(e-e a )(e-e v ) = J 2 



(10) 



Since we consider a low-energy spectrum |e| <C A u , the tra- 
jectories should pass close to the vortex center for the spec- 
trum modification (TlOb to be effective. Then, the interaction 



of surface and vortex states is determined by the overlap in- 
tegral J sa Aexp(— a/£), where a = a/ cos9 p and a is the 
distance from the vortex to the surface. Taking a certain point 
at the surface (see point A in Fig|2^) of the superconductor 
one can obtain a relation between the angles and impact pa- 
rameters of trajectories passing trough this point as follows 
b = d s'm(9 — 9p). Thus the energy of vortex core states can be 
written as e v = M(a/£)Arj sm(9 — 9 P ). Then, from Eq.dTOt 
we obtain the spectrum transformation shown qualitatively in 
Fig|2]for the particular case of 9 = 0. It is easy to see that for 
equal values of vorticity and chirality (Fig|2j)) there appears 
a minigap in quasiparticle spectrum near the Fermi level and 
therefore the zero-energy DOS is suppressed. On the other 
hand, in the case of opposite vorticity and chirality (FiglJj;) 
there is no minigap and the DOS is not suppressed. 



(a) 



A J 




incident QP 
M -3j 



vortex 





surface surface ^ 
vortex 



FIG. 2: Sketch of QP trajectories forming surface and vortex states 
(a) and qualitative plot of spectrum transformation due to the inter- 
action of surface and vortex states in the case of a chiral p-wave su- 
perconductor; vorticity and chirality have (b) equal and (c) opposite 
values. 

In a d+is- and d+id- wave superconductor, the interaction 
between vortex and surface states can also lead to noticeable 
effects, which will be discussed later in the present paper. 

Recently, it was pointed out that tunneling of quasiparti- 
cles into vortex core states leads to a resonant enhancement 
of subgap conductance of normal metal/superconductor (N/S) 
junction 40 . In the case of chiral superconductors, such a tun- 
neling effect can lead to either stimulation or suppression of 
conductance, depending on the direction of vorticity. We will 
show that if vortices are located far from the N/S interface, the 
conductance follows the behaviour expected from the Doppler 
shift approach. On the other hand, when the distance from 
the vortex to the interface becomes comparable with coher- 
ence length £ the tunneling into vortex core states comes into 
play, leading to the peculiar nonmonotonic conductance de- 
pendence on the vortex coordinate with respect to the super- 
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conducting surface. 

This paper is organized as follows. In Sec. [II] we give an 
overview of the theoretical framework which is employed in 
this work, namely a Bogolioubov approach and a quasiclas- 
sical Eilenberger approach. In Sec. [TTTJ we present our main 
results for the influence of magnetic field on bound surface 
states spectra and local density of states near the surface. We 
discuss the transformation of surface states in the Miessner 
state of superconductor as well as the effects of interplay be- 
tween surface and vortex core states. We give our conclusions 
in Sec. [TV] 

II. THEORETICAL APPROACH 

Our further considerations are based on the Bogoliubov-de 
Gennes (BdG) equations for particle- (u) and hole-like (v) 
parts of the wave function, which have the following form: 

1 / e \ ^ 

- - — (p A u + Av = (s + e F )u, 

2m V c / 




p+-A) v + AJu = (e-e F )v. (11) 



Here A is a gap operator, A is a vector potential, p — 
—i(d/d x , d/dy), and r = (x, y) is a radius vector in the plane 
perpendicular to the anisotropy plane. Hereafter we assume 
the Fermi surface to be cylindrical along the z- axis and con- 
sider a motion of QPs only in xy plane. 

In case of unconventional superconductors, the gap po- 
tential A is a non-local operator, so the BdG system effec- 
tively becomes a very complicated integro-differential equa- 
tion. Another complexity arises from the broken spatial in- 
variance of the superconducting gap in presence of vortices 
near the N/S interface. A simplification can be obtained if 
one considers a quasiclassical approximation, assuming that 
the wavelength of quasiparticles is much smaller than the su- 
perconducting coherence length (see e.g. Refill). Within 
such an approximation, QP move along linear trajectories, 
i.e. straight lines along the direction of QP momentum 
n = kFfc^ 1 = (cos Op, sin 9 P ). Generally, the quasiclassi- 
cal form of the wave function can be constructed as follows: 
(u, v) — e lkF " r (t/, V), where (U(r), V(r)) is a slowly vary- 
ing envelope function. Then the system (fTTT i reduces to a 
system of first-order differential equations along the linear 
trajectories defined by the direction of the QP momentum 
n = kpkp 1 — (cos 9 p , sin 9 p ). Introducing the coordinate 
along trajectory x' = (n • r) = r cos(9 p — 9) we arrive at the 
following form of the quasiclassical equations: 

{-ihv F d x > + v F ■ ^A) U + AV = eU 
(ihv F d x , + v F --A)V + A^U = eV, (12) 

where the Fermi velocity is vf = nhk F /m. The pairing 
potential in Eq. ( TTZb may generally be written as 

A(r,0 p )=A(0 p )¥(r), (13) 



where A(9 p ) describes the orbital symmetry of the super- 
conducting order parameter in momentum space, while (r) 
describes its spatial dependence both magnitude- and phase- 
wise. 

The local DOS (LDOS) can be expressed through the eigen- 
functions of the BdG equation (fTTT i in the following form 32 : 

N(e,r) =^K(r)| 2 5(e-e„), (14) 

n 

where it„(r) is electron component of quasiparticle eigen 
function corresponding to an energy level e„. The eigenfunc- 
tion has to be normalized: JJ^° \u n (r)\ 2 + \v n (r)\ 2 d 2 r = 1. 

We will also later employ a quasiclassical Eilenberger ap- 
proach to study the spatially resolved DOS. Let us here sketch 
the framework of the treatment which makes use of the Eilen- 
berger equation, following the notation of Refs . 35 ' 36 . It is now 
convenient to solve the Eilenberger equation along trajecto- 
ries along the Fermi momentum, and introducing a Ricatti- 
parametrization for the Green's function 36 . In this way, one 
obtains 2 ^ 

hv F d x >a{x') + [2u) n + A f a(x')]a(x') - A = 0, 
Hv F d x ,b(x') - [2Q n + Ab(x')]b(x') + A f = 0, (15) 

where ia>„ = iui n + mv F • v s is a Doppler-shifted Matsubara 
frequency and 

v s = J-(W*--A > ) 

2m \ c J 

is a gauge-invariant superfluid velocity where $(r) is a gap 
function phase: ^(r) = l^le 1 *. The LDOS may be expressed 
through the scalar coherence functions a and b as follows 20 

J 2tt ( 1 + ab J iuj n ^e+id 

where e is the quasiparticle energy measured from Fermi level 
and S is a scattering parameter which accounts for inelastic 
scattering. 

To investigate the transport properties of N/S junction, we 
employ an approach similar to what was used in work by 
Bardeen, Tinkham and Klapwijk 33 . The expression for the di- 
mensionless zero- bias conductance of the N/S junction mea- 
sured in terms of the conductance quantum e 2 /(7rS) can be 
written as follows: 

G=%/ / [I -R n (9 ) + R a (0 )] cos Q d9 , (17) 

Z J-tt/2 

where R n (9o) and R a (0o) are the probabilities of normal 
and Andreev reflection respectively, 9q is the incident angle: 
kp = k F (cos 9o, sin6*o), characterizing the propagation di- 
rection of quasiparticles, coming from the normal metal re- 
gion. The Sharvin conductance G s h = k F L y /ir equals the 
total number of propagating modes determined by the chan- 
nel width L v . 
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The problem of quasiparticle scattering at the N/S interface 
is formulated within the BdG theory dTTb . An interfacial bar- 
rier separating the N and S regions can be modeled by re- 
pulsive delta function potential W(x) = Wq6{x), parameter- 
ized by a dimensionless barrier strength Z = Wq/Hvf- The 
boundary conditions at the N/S interface then read: 

[/(0)]=0, [d x f(Q)] = (2k F Z)f(0), (18) 

where / = (u, v) and [f(x)] = f(x + 0) - f(x - 0). 

Considering a zero-bias problem we will have to ana- 
lyze only zero-energy excitations with e = 0. For wave 
functions in S region corresponding to subgap quasiparti- 
cles, the following representation can be used: (U, V) — 
e C ( e i(v+*)/2 ie -i(n+*)/2y W here C = CM) and r) = 
r](s, b) are real-valued functions. Then, the quasiclassical 
equation Sl2\ can be written as follows: 

d x ,r] + 2|A| cost? + 2e d = 0, (19) 
9 a /C + 2|A|sinr? = 0. 

where £d(r) = Mcfv s is a Doppler shift energy. For 
wave functions decaying at the different ends of trajectory 
(U, V)(x' = ±oo) = from Eqs.fO we obtain: 

r)(x' = ±oo) = ±tt/2. (20) 

The boundary conditions ( TT8l model the specularly reflect- 
ing N/S interface, coupling the waves with wave vectors kp = 
fcF(cos£?0;Sin(9o), an d kV = fci?[cos(7r — #o),sin(7r — 9q)}. 
Therefore if the incident electron wave is Ui = e lkFr , then the 
reflected electron u r and hole v r waves will have the form 

u r = U r e lk '* r , v r = K.e lkFr , 

where U r and V r are the envelope functions. Thus, each point 
(0, y) at the N/S interface lies on the intersection of two qua- 
siclassical trajectories, characterized by the angles 9 p = 9q 
and dp = 7r — 6>o- Let us denote the distribution of phases 
r){x') along these trajectories as ij + (x / ) and r]-(x') corre- 
spondingly. Using the boundary conditions we obtain the fol- 
lowing expression for the conductance^: 

G=-£ / g(y,0 o )cos9 o d0 o dy, (21) 

1 J -Ly/2 J -ir/2 

where g(y, 8 ) is given by 

g(y,0o) = ^ = — , (22) 

v ; (Z 4 + Z 2 )|l- e v| 2 + l 

with Z = Z/ cos 9q and p(y, 9o) = r;_ — r/ + is determined by 
a difference of phases f]-(x') and r) + (x') at the intersection 
point (0, y). To evaluate the conductance, one needs to find 
the factor e ip in Eq.(l22ii and then the reflection probabilities 
by solving numerically Eq.([T9l with the boundary conditions 
in Eq. (gDJ. 



UI. RESULTS 

To illustrate the basic effect of how the interplay between 
the Doppler-shift and the time-reversal symmetry breaking of 
the superconducting order parameter is manifested, we con- 
sider a situation where an external magnetic field is applied 
near the surface of the superconductor along the z-axis, thus 
inducing a vector potential A in the superconductor which 
drives the shielding supercurrent. In order to proceed analyti- 
cally, we make the simplifying assumption that the superfluid 
velocity field is nearly homogeneous and that the spatial varia- 
tion of the superconducting order parameter near the interface 
is small. Choosing a real gauge, we then find that the Ricatti- 
functions a and b in Eq. (T5[ may be written asj22i22 

a{9) = s(6)A(6), b(9) = s(6)A*(9), 

8(6) = l/[w n (0) + VK(0)) 2 + |A(0)| 2 ], (23) 

where ui n depends on 9 through the Doppler-shift. To evaluate 
the LDOS in Eq. (TToT l at the surface, we need to take into 
account proper boundary conditions at x = 0. Assuming an 
impenetrable surface with perfect reflection, these boundary 
conditions read: 

a surf ace(60 = a(TT - 9), b sm{ace (9) = b{9). (24) 

Inserting this into the expression for the LDOS, we obtain 

Af(e) = 2Ref(- , 1 -1. (25) 

w l\l + a(n - 9)b(9) / ) iu n ^e+i8 

The (...) denotes angular averaging, which we restrict to an- 
gles —ir/2 < 9 < tt/2 due to the surface. It may be shown 
that for a chiral p-wave superconductor 22 , the zero-energy 
DOS at the surface reads 

AT(o) = i + -^1. + . . . , (26) 

while for pure s- or rf-wave superconductors one finds 

AT(0) = d + C 2 v% + ..., (27) 

where C\ and C-2 are arbitrary constants. From numerical in- 
vestigations of Eq. d25b at e = 0, we find that the zero-energy 
DOS may quite generally be written as 

AT(0) = C x + C 2 v sy + ..., (28) 

whenever the superconducting order parameters i) break time- 
reversal symmetry and ii) support the presence of subgap 
surface-bound states. This is the case both for the p x +ip y - 
wave pairing which is believed to be realized in Sr2Ru04, as 
well as the c?+is-wave and d + id-wave pairings that are rel- 
evant for the cuprates. In particular, tunneling spectroscopy 
measurements have indicated the presence of such a time- 
reversal symmetry breaking order parameter near surfaces by 
a split zero-bias conductance peak that was observed in the 
absence of an external field in several experiments^. 

In Ref. |35l it was pointed out that the neglect of the gradient 
term in the Eilenberger equation is expected to be a reasonable 
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approximation as long as the Doppler-shift energy mvp ■ v s 
is small compared to the local gap energy A(9). This approx- 
imation would then fail close to the vortex core or gap nodes 
of A(0). Nevertheless, in the model case of spatially homoge- 
neous gap function and superfluid velocity field the gradient 
terms in Eilenberger equation can be neglected in the whole 
range of Doppler shift energies. However considering a model 
situation the above discussion nevertheless serves to illustrate 
our main qualitative argument: namely, that chirality-sensitive 
effects should be expected in superconductors with order pa- 
rameters that i) break time-reversal symmetry and ii) support 
the presence of subgap surface-bound states. We now proceed 
to discuss the cases of p x + ip y -wave and d + is(d)-wave pair- 
ing in more detail, since these are relevant to actual materials. 



A. Surface states in p + ip and d + is(d) superconductors 
under the influence of magnetic field. 

In Figf3] we show numerical plots of the surface LDOS 
at the Fermi level given by Eq.d25l for the chiral p-wave 
(Fig|3k), d + is-wave (FigfJJ)) and d + id-wave (Figj3j;) cases 
in a wide domain of superfluid velocities. The structure of gap 
functions is chosen in the form of Eqs. (H}-© and the param- 
eter characterizing inelastic scattering in Eq.(fT6b is chosen as 
5 = O.lAo. We introduce the following notations for the dif- 
ferent critical velocities: v c — Ao/(hkp), v cs = \ A s \/(hkp) 
and v cd = \A d \/(hk F ). 

As seen, the surface LDOS has sharp peaks at a certain 
value of the superfluid velocity in all cases. We will show 
below that peaked structure of LDOS is provided by bound 
surface states. Another contribution to the LDOS comes from 
the delocalized states corresponding to the continuous part of 
QP spectrum. A delocalized state with zero energy e = 
exists provided (i) \v sy \ > v c in case of chiral p-wave super- 
conductor, (ii) \v sy \ > v cs and (iii) \v sy \ > v c d in case of 
d + is-wave and d + id-wave superconductors correspond- 
ingly. Condition (i) is unlikely to be realized because it means 
that the superfluid velocity is larger than the critical depair- 
ing value. Conditions (ii) and (iii) can be realized, because 
the values v cs and v c d can be well below the critical depairing 
velocity if the amplitude of additional order parameter com- 
ponents is small enough. 

To analyze the contribution to LDOS provided by the bound 
surface states we will consider the domain of low energies 
|e| <C Ao- By neglecting small deviations of the electron and 
hole momentum, the normalized wave function of QP local- 
ized near the boundary can be written as 








£1 cos6> p | 



where (k x ,k y ) = kF(cos9 p ,s'md p ). This wave function de- 
cay in the superconducting side x > at a characteristic lo- 
calization scale £ is given by £ = hvp/Ao for chiral p-wave 
and £ = hvp/\Ao sin(2# p )| for d + is- and d + id- wave 
superconductor correspondingly with gap functions given by 
Eqs. (l2l3b . The spectrum of the Andreev bound states, shifted 
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FIG. 3: Plot of the normalized zero-energy LDOS N(0) for (a) p- 
wave superconductor with x — L (b) d+is-wave case with A s = 
0.4Ao and (c) d+id-wave case with Ad = 0.4Aq. Dashed lines are 
guides for eyes: the vertical ones denote positions of LDOS peaks 
and horizontal ones correspond to the level of normal metal DOS 
No. 



by the superfluid velocity, is given by 

e a = x^-aky/kp + hv sy k y (29) 
for the p-wave case, 

e a = A s sgn(fc v ) + hv sy k y (30) 

for d + is-wave case and 

e a = A d sgn(fc y ) cos(26» p ) + hv sy k y (31) 

for d + id- wave case. Consequently, the contribution from 
Andreev bound states to the zero-energy LDOS at the surface 
of a chiral p-wave superconductor is given by 



N a = N - 



1 



Vsy/Vc + xY 

where N = m/(27rfi. 2 ) is the normal metal LDOS per one 
spin direction. For a chiral d + is superconductor, the be- 
haviour of the LDOS is more complicated. Assuming that 
A s > 0, we obtain that the LDOS is zero for v sy > 
—As/hkp- Otherwise, it is given by 

v cs v c 



N a = 4/Vq- 
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On the contrary, for the d + id case the LDOS is zero if v sy > 
Aa/ (hkp) (for A,^ > 0) and otherwise it is given by 

n "= n °£i ( i+ i vcd V 

It can be seen that these contributions to LDOS have peaks 
at v sy = v c for p-wave case. For d + is-wave and d + id- 
wave superconductors the peaks are positioned at v sy = 
— sgn(A s )v cs and v sy = sgn(Ad)v c d correspondingly. Even 
though the position of the peaks are different, the dependen- 
cies of the surface LDOS on the superfiuid velocity (and con- 
sequently on magnetic field) are very similar for d + is and 
d + id- wave superconductors. Therefore, it might be difficult 
to distinguish which case is realized experimentally. 

On the other hand the considered model with a spatially ho- 
mogeneous gap function ^(r) = 1 is adequate only when the 
applied magnetic field is not too large. When the magnetic 
field is large enough, it breaks the Meissner state and gener- 
ates vortices near the surface of superconductor. Therefore, 
we investigate the influence of vortices on the LDOS distri- 
bution near a superconducting surface as well as on the con- 
ductance of normal metal/superconducting junctions. We will 
show that vortices have different effect on the conductance in 
d + is and d + id cases. 



B. Interplay of vortex and surface states in chiral 
superconductors. 

A chirality sensitive LDOS transformation due to vortices 
situated near the surface of a chiral p-wave superconductor 
was considered in Ref|H It was shown that depending on the 
chirality and vorticity value, the surface LDOS near is either 
enhanced or suppressed upon decreasing the distance from the 
vortex to the surface. In case of d + is(d) superconductors the 
transformation of LDOS profile is also sensitive to the value of 
vorticity. Similar behaviour is expected for a conductance of 
normal metal / chiral superconductor junction in the presence 
of vortices. 

To investigate the influence of a single vortex on the LDOS 
profile and conductance, we assume that at x > (supercon- 
ducting region) the coordinate dependence of the order param- 
eter may be written as follows: 

tf (r) = e 1 *. (32) 

Here, we consider a model situation where the magnitude of 
the order parameter is constant. The phase distribution <£>(r) 
consists of a singular part & v (r) = arg(r — r v ) and a reg- 
ular part < I > r(r), determined by the particular metastable vor- 
tex lattice configuration realizing near the boundary. We as- 
sume that the regular part of the phase distribution is $ r (r) = 
—arg(r — r av ) corresponding to the image vortex situated at 
the point r av = (—2a, 0, 0) behind the N/S interface. 



1. p + ip wave 

In Fig 0] we show the LDOS profile near the surface of a 
chiral p-wave superconductor in the presence of a single vor- 
tex, positioned at some distance a from the surface. When the 
vortex is positioned far from the surface a > 2£ the LDOS 
profile follows the behaviour, expected from the picture of lo- 
cal Doppler shift 22 . Depending on the relative value of vortic- 
ity and chirality, the surface LDOS is either increased (Figj4^) 
or decreased (Figj4j)). An analytical estimate with the help of 
spectrum Eq. (l33l l yields a following estimation of the am- 
plitude of LDOS peak in Figgji: AN/N = (1 + Mxa)' 1 . 
At smaller distances a < 2£, the behaviour of LDOS changes 
drastically. In the case of opposite vorticity and chirality, the 
surface LDOS grows at a < 2£, obviously due to the over- 
lapping with the peak of vortex core states. In case of equal 
vorticity and chirality the same overlapping occurs, but on the 
contrary it leads to reduction of DOS, as it was discussed in 
Sec. ID The peak of the LDOS at the surface discussed in 
Refl22l transforms into a dip-and-peak structure as the vortex 
comes close to the surface. 




FIG. 4: (Color online) Plot of the normalized zero-energy LDOS 
N(0) in the presence of a vortex near the surface of a chiral p-wave 
superconductor, (a) and (c) correspond to equal vorticity and chiral- 
ity, (b) and (d) correspond to opposite vorticity and chirality. The 
distance from vortex to the surface is a — 2£ for (a) and (b) and 
a = £ for (c) and (d). 

This is illustrated in FigOa), where we plot the LDOS at 
the surface point (0, 0), which is the nearest point to the vortex 
in FigH] At large distances a ^> £ the LDOS is a monotonic 
function of a, either increasing or decreasing depending on 
the relation between vorticity and chirality. At smaller dis- 
tances a < 2£, the extremum of LDOS appears. In the case of 
opposite vorticity and chirality [lower curve in FigfSJa)], the 
surface LDOS grows at a < 2£, due to the overlapping with 
the peak of vortex core states. In case of equal vorticity and 
chirality [upper curve in FigfSJa)] the same overlapping leads 
to reduction of LDOS. 

To investigate the influence of vortices on the transport 
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properties of normal metal/ chiral p- wave superconductor 
junction we solve the generic problem of the influence of a 
single vortex near the N/S surface on the zero-bias conduc- 
tance of the junction. A numerical plot of the conductance 
G as a function of a distance of vortex to the junction inter- 
face is shown by the solid lines in Figj5] (b) for equal (upper 
curve) and opposite (lower curve) values of chirality and vor- 
ticity. The conductance is normalized to the value of Sharvin 
conductance G s h — kpL y /-K. 



(b) „, G/G, t 
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FIG. 5: (a) Plot of the normalized zero-energy LDOS at the point 
on the surface which is closest to the vortex core. Different curves 
correspond to different vorticities. (b) Plot of the vortex-induced con- 
ductance in chiral p+ip-wave superconductor for equal and opposite 
values of vorticity and chirality. The strength of interface barrier is 
Z — 5. Large distances asymptotics for N and G are shown by dash 
lines. 

At large distances a ^> £ an analytical estimation of con- 
ductance can be obtained by using a local Doppler shift ap- 
proximation on the quasiparticle spectrum. Indeed, the modi- 
fication of the surface states energy due to a supercurrent flow- 
ing along the boundary of superconductor can be written as 



(xA 



hV S yk F ) 



ky / kf 



(33) 



where k y is a quasiparticle momentum along the surface, 
X = ±1 is a chirality value and v sy — (MH/m)a/(y 2 + a 2 ) is 
a projection on the surface plane of superfluid velocity gener- 
ated by the vortex and image antivortex, M is vorticity value 
and m is the electron mass. It follows from Eq.d33l that the 
Doppler shift effect leads to a change in the slope of anoma- 
lous branch. It is easy to obtain that in this case the function 
g(y, 6q) in expression Eq. ( f22T > takes the following form: 



A(Z^ + Z 2 )(e a /A Q ) 2 + l 



(34) 



The straightforward integration in Eq.(t2Tb yields G = Gq + 
SG, where Go = G s h{it/Z 2 ) is the conductance without vor- 
tex and 



5G/G sh = ± 



2tt£ 
Z 2 L V 



arctan(L !; /2a) 



(35) 



is a vortex-induced conductance shift, where the upper (lower) 
sign corresponds to equal (opposite) vorticity and chirality. 

At distances smaller than 2£, an extremum of the conduc- 
tance appears. Upon placing the vortex closer to the surface, 



an opposite effect occurs: one obtains a conductance suppres- 
sion instead of enhancement and vice versa. The origin of the 
conductance extremum is a tunneling of quasiparticles into the 
vortex core states, or in other words, the overlapping of vortex 
and surface bound states. Comparing FigsOa) and [2b) one 
can see that the conductance in general follows the behaviour 
of the surface DOS. 



2. d + is and d + id wave. 

In chiral d+is and d + id superconductors the LDOS trans- 
formation appears to also be vorticity sensitive. In the Fig(6] 




FIG. 6: (Color online) Plot of the normalized zero-energy LDOS pro- 
file N(0) in the presence of vortex near the surface of chiral d + is- 
wave superconductor with A s = 0.2Ao. (a),(c) and (b),(d) corre- 
spond to different vortex orientations with respect to z axis. The 
distance from vortex to the surface is a — 4£ for (a) and (b) and 
a = 2£ for (c) and (d). 

we show the profile of zero-energy LDOS in the case when 
the vortex is placed at a distance of a — 2£ from a flat bound- 
ary of a d + is-wave superconductor characterized by a gap 
function in momentum space given by Eq.©. In this section, 
we use the notation £ = hvp /Aq. 

One can see that for one sign of vorticity the surface LDOS 
shows two peaks which are symmetric with respect to the vor- 
tex position. As we have shown above, the large peaks in 
surface LDOS appear when the energy coincides with the po- 
sition of bound state level. For a different sign of the vorticity, 
there are no surface states at the Fermi level and the LDOS 
along the surface is a flat function. A non-zero level of LDOS 
in this case is provided by inelastic scattering which leads to 
the smearing the QP energy levels. Applying a local Doppler 
shift approach, which holds if the distance from vortex to sur- 
face is rather large (a ^S> one can interpret the results shown 
in Fig.©. 

The coordinates y* of surface LDOS peaks can be es- 
timated from the relation v sy = A s /pp, where v sy = 
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FIG. 7: (Color online) Plot of the normalized zero-energy LDOS pro- 
file N(0) in the presence of vortex near the surface of chiral d + id- 
wave superconductor with Ad = 0.2Ao (a),(c) and (b),(d) corre- 
spond to different vortex orientations with respect to z axis. The 
distance from vortex to the surface is a = 4£ for (a) and (b) and 
a = 2£ for (c) and (d). 



(Mh/m)a/(y 2 + a 2 ) is a projection on the surface plane 
of superfluid velocity generated by the vortex with vortic- 
ity M and image antivortex. It can be seen that for a > 
£(A /|A S |) the peak is situated at y* = i.e. at the 
surface point nearest to the vortex. Otherwise, we obtain 
y* = ±ayl — (a/£)(|A s |/Ao). Comparing this estima- 
tion with the numerical results in Figj6] one observes a mi- 
nor difference. For example, it follows from the estimation 
that the LDOS peaks should be positioned at y* — ±2.4£ for 
A s = 0.2A , but in Fig|6]they are located at y* — ±2.0£. 
This discrepancy can be attributed to the complex shift of the 
energy e — > e + iS due to the effective scattering parameter 
5 = O.lAo which was used in the numerical calculations. If 
we increase the distance from vortex to surface a, the LDOS 
peaks will merge when a > £(Ao/|A s |) (see Fig|6]l. 

In FigfT] we show the LDOS profile modulated by a vor- 
tex placed at a distance of a = 2£ from a flat boundary 
of d + id superconductor. The structure of the gap func- 
tion was chosen in the form ©. Applying the approach 
based on the local Doppler shift we obtain the similar ex- 
pression for the coordinates of the peaks of surface LDOS: 
y* = ±ayl — (a/£)(|Ad|/Ao). For the particular values of 
parameters A^ = 0.2Ao and a = 2£ this estimation yields 
y* = ±2.4£, which is much less than obtained from numer- 
ical plot in Figj7] y* ±4£. This discrepancy can also be 
attributed to the effect of inelastic scattering, which appears 
to have a larger effect in d + id- wave case than in discussed 
above d + is-wave case. 

A numerical plot of the N/S junction conductance as a func- 
tion of distance from vortex to surface is shown in Figj8]for the 
d + is and d + id cases. The conductance is normalized to the 
Sharvin conductance G s h = kpL y /n. Comparing Fig|8^ and 
FigJHJ) one can see that the conductance behaviour is qualita- 



tively different for s and d- wave symmetry of the additional 
gap function component. For d + is- wave, the conductance 
has a sharp peak for one vortex orientation (upper curve in 
FigEh) and it is a flat function of a for another vortex orien- 
tation (lower curve in Fig|8^). The origin of the conductance 
enhancement is a formation of Andreev bound states at the 
Fermi level which are localized near the superconducting sur- 
face. As was discussed in the Introduction (see the Fig. []J>), 
the zero-energy Andreev bound states can appear only for a 
certain direction of superfluid velocity flowing along the su- 
perconducting surface and if the value of the superfluid veloc- 
ity is larger than a critical value \v sy \ > \A S \/ (Hkp). For a 
high interface barrier Z> 1 applying an approximate analyt- 
ical expression (l35l l we find that a sharp increase of conduc- 
tance in FigJSh can be described by the following expression 



G/G s h — 



16tt £ A 



a \3/2 
a* 



where a* = £(A /|A S |) and A ~ 1. Otherwise, if a > a* the 
conductance is much smaller since Z>1: 



G/G s h 



Ao 



4 



When the distance a is decreased further, the conductance is 
suppressed (see Fig[8^, upper curve). The decrease of con- 
ductance can be attributed to the gap at the Fermi level which 
appears due to the interaction of vortex and surface states in 
a similar way as for the p + ip- wave case discussed in the 
previous section. 

In a d + id-wave superconductor, zero-energy Andreev 
bound states may exist even in the absence of vortex. An 
asymptotic value of the conductance Go at a ^ ^ can be ob- 
tained using the expression ( f35l > as follows 



Gq/ G s 



Ad\J 2V2Z 2 ' 



When the vortex approaches the superconducting surface, the 
conductance is either suppressed (lower curve in FigJHJ?) or 
slightly enhanced (upper curve in FiglS}}). This behaviour can 
be understood by again using the Eqs.(|35ll with the Doppler 
shifted spectrum of Andreev bound states ([6j. The decrease 
(increase) of conductance corresponds to the transformation 
of spectrum shown qualitatively in FigQJ by dash (dash- 
dotted) lines. It is possible to obtain an analytical expression 
for the vortex-induced conductance shift at a ^> £ in the fol- 
lowing form: 



SG/G S 



± 



2Z 2 



A d 



arctan 



where the upper and lover signs correspond to the different 
vortex orientations. As the vortex approaches the surface fur- 
ther, there appears an extremum of the conductance. Such be- 
haviour can be explained by a conductance enhancement due 
to the tunneling of QP into the vortex core states, discussed in 
Refkoi A sharp decrease of the upper curve in Fig [8}? can be 
attributed to the opening of an energy gap at the Fermi level 
due to the interaction of vortex and surface states. 
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FIG. 8: Plots of the vortex-induced conductance in case of (a) chiral 
d+is superconductor for A s = 0.2Ao and (b) d+id superconductor 
for Ad = 0.2Ao. The strength of interface barrier is Z — 5. Differ- 
ent curves on each plots correspond to different vortex orientations 
with respect to the z- axis. 



IV. SUMMARY 

In summary, we have investigated how the tunneling con- 
ductance and the local density of states (LDOS) in supercon- 
ductors are affected by the influence of an external magnetic 
field when the superconducting order parameter (OP) breaks 
time-reversal symmetry (TRS). This is directly relevant for 
both Sr2RuC>4, where chiral p + ip-wave pairing is believed to 
be realized, and for the high-T c cuprates, where a d + is- or 
d + id-wave OP has been suggested to exist near surfaces. 
In addition to breaking TRS, all of these OPs feature sur- 
face bound zero-energy states at surfaces under appropriate 
circumstances (e.g. a dominant c?-wave OP in the d + is-wave 
case). 

We have shown how the Doppler-shift conspires with an 
interaction of vortex and surface states to produce a consid- 



erable qualitative modification of both the tunneling conduc- 
tance and the LDOS. When the vortex is located at distances 
well above a coherence length £ from the surface, the Doppler- 
shift produces an enhancement or suppression of the LDOS 
depending on the relative sign of the vorticity and the chirality 
of the superconducting OP. This effect may be directly probed 
by first applying an external magnetic field in a direction while 
measuring the LDOS and then reversing the field direction and 
measuring again. When the vortex is located very close to the 
surface (a distance of order £ or smaller), there is an overlap 
between the vortex and surface states which effectively cause 
a dramatic change in the tunneling conductance and LDOS. 
This effect is also sensitive to the relative sign of the vortic- 
ity and the chirality of the superconducting OP. The overlap 
between these two sets of states results in either a strongly 
enhanced or suppressed tunneling conductance/LDOS at zero 
bias voltage/zero energy. 

We have demonstrated the aforementioned effects both 
qualitatively and quantitatively for p+ip-, d + is-, and d + id- 
wave symmetries. Experimentally, the distance from the sur- 
face to the closest vortex can be altered by modifying the field 
strength. All of our predictions should be possible to test ex- 
perimentally with present-day techniques. 
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